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Abstract 

The quantum state of a particle can be completely specified by a po- 
sition at one instant of time. This implies a lack of information, hence 
a symmetry, as to where the particle will move. We here study the con- 
sequences for free particles of spin and spin 1/2. On a cubic spatial 
lattice a hopping equation is derived, and the continuum limit taken. Spin 
leads to the Schrodinger equation, and spin 1/2 to the Weyl equation. 
Both Hamiltonians are hermitian automatically, if time-reversal symmetry 
is assumed. Hopping amplitudes with a "slight" inhomogeneity lead to the 
Weyl equation in a metric-affine space-time. 



For a state vector which is a superposition of eigenvectors of an observable, 
the probability for measuring an eigenvalue is \ipi\ 2 if ipi is the corresponding 
amplitude. The prevailing view about this rule is that it must be taken as 
an independent Statistical Axiom. Several authors, however, have recently 
presented arguments to show that the relation of amplitudes and probabilities 
is contained already in the non-statistical axioms of quantum mechanics 
H, ||, U, thus raising an issue of potential interest relating to the difficulties in 
defining a satisfactory scalar product for state vectors in General Relativity 
||. Clearly, there were assumptions to be made. In Jl|, |2|, |3| the normalization 
of state vectors in the usual sense was anticipated; in || it was assumed that 
in an isolated two-state system there are no "decaying" modes of evolution. 
The present paper intends to weaken the latter assumption (thus generalizing 
the conclusion) from "non-decaying" to "non- singular" . Instead of two-state 
systems, we here consider free particles of spin and of spin 1/2. For these, 
the unitarity of time evolution can be inferred from merely 

• the superposition principle 

• symmetry under translation, rotation, and time-inversion 

• the "logic" of position eigenstates. 

The argument is very simple for a scalar particle living on the sites of a 
cubic lattice. The particle's position at the instant of preparation already 
specifies the state vector (up tp a phase), so there is an unavoidable inde- 
terminacy as to where the particle will go next. On a lattice, all nearest 
neighbours are equivalent destinations of the hopping process. This amounts 
to a symmetry which is not present in classical motion. Indeed, only by 
forming superpositions — the unique quantum-mechanical concept relating to 
undecided alternatives — can a particle "move" out of a position eigenstate. 
The symmetry among the nearest neighbours in this process allows to re- 
construct the one-particle Hamiltonian without any reference to the classical 
equations of motion, as elaborated already in 0. In the present paper, more- 
over, we avoid anticipation of a scalar product, invoking higher symmetries 
and general superpositions instead. 

The argument takes its most suggestive form in the Heisenberg picture. 
Let a be the spacing of a cubic lattice and let 

\x,t) x = an, n integer 



2 



be an eigenstate of the position operator, prepared at the time t. When this 
same state vector is expressed in terms of position eigenstates referring to 
a time infinitesimally different, t' — t + dt, we expect to be dealing with a 
superposition of nearest-neighbour states, i.e. 

\x, t) = a\x, t + dt) + V] P\x + ah, t + dt) (1) 

nearest 
neighbours 

where a and (3 are some complex numbers dependent on the time step dt. 
For dt — ► we expect a — > 1 and (3 — > 0; hence, 

a = 1 + edt + C(dt 2 ) /3 = K dt + 0{dt 2 ) 

Thus (HD amounts to the differential equation 

— — \x, t) = e\x, t) + ^2 n\x + ah,t) (2) 

^ nearest 

neighbours 

Invoking the superposition principle again, we now consider a general state 
vector \ip) in the Heisenberg picture. It is time-independent, or rather "time- 
less" [0 as a matter of concept, but \i[)) can be expressed in a basis of position 
eigenstates referring to any time t' . Thus 

lattice 
sites 

where ip{x,t') denote some complex coefficients. Differentiating (|3|) with 
respect to t', using the linear independence of \x, t') at any fixed t', and using 
equation (0), we conclude that 

^ijj(x,t') = eip(x,t') + J2 ki/j(x + ah,t') (4) 

nearest 
neighbours 

In order to take the continuum limit, which physically could mean something 
like a — > api anc k rather than a — > 0, we expand 

+ an, t') = + ahiViil){x, t') + la 2 h i h j V i V j 'ip(x, t') + C(a 3 ) 

In the sum over nearest neighbours in equation (Q) the first-derivative terms 
cancel out while second-derivative terms remain to form the Laplacian of i/j. 
Thus, changing t' into t again, 

d 

—ip{x, t) = (e + 6K)ip(x, t) + a 2 KAip(x, t) + 0(a 3 ) 
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When terms of 0(a 3 ) are omitted, the general solution of the equation is 

^(f ; t) = e {e+6K)t J ^($e^V KoV 'd 3 p (5) 

where ip(p) is the Fourier transform of the wave function at t = 0. Assuming 
that quantum mechanics has a time-reversal symmetry we can see from ([5]) 
that hopping parameter k must be purely imaginary. This is because ^tn ^ 
would imply an evolution factor exponentially increasing like 

exp (jKftla 2 ]? 2 ^!) 

in either the positive or negative time direction. Thus any wave function at 
t = with a fall-off in momentum space slower than Gaussian would run into 
a singularity immediately at either t < or t > 0. The other evolution factor 
in equation (|5|), e ( £+6K ^, is independent of the one-particle state vector \?p). It 
would, however, multiply to the nth power in a state with n particles. Hence 
it would lead to a singularity similar to the above for any superposition of 
n-particle states with a fall-off in n slower than exponential. Therefore e, 
too, must be purely imaginary. Thus e (e+6K - )t can be absorbed in a universal 
redefinition of the phase of the wave functions. 

The argument does not depend on the number of spatial dimensions. It 
can be used to corroborate the "no decay" assumption of 0, relating to the 
propagation of a quantum particle through a "channel" connecting "cavities" . 
If the "channel" is modelled by a one-dimensional array of lattice sites, then 
in the continuum limit the propagation is unitary automatically. 

Let us now consider a spin 1/2 particle living on the sites of a cubic 
spatial lattice. A similar scenario was investigated by Bialynicki-Birula || in 
connection with cellular automata; our approach differs by letting time run 
continuously, and by not anticipating the unitarity of the time evolution. Let 

\x,s,t) x = an s = ±| 

be an eigenstate of position and spin prepared at the time t. In terms of 
state vectors prepared at f = t + dt we expect a superposition involving 
nearest-neighbours and spin-flips. As in the scalar case, let nearest-neighbour 
directions be given by unit vectors h with n = ±1, ±2, ±3, where — n = —n. 
For convenience, let = denote the unshifted position. We then expect a 
hopping equation of the from 

3 

\x, s,t) = c nss f\x + an, s', t + dt) (6) 

n=— 3 
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where we adopt the summation convention with respect to the spin. In the 
limit dt — > the RHS must tend to the lhs, so 

c 0ss > = 5 SS > + H 0ss > dt + G(dt 2 ) 

W = H nss ,dt + (D(dt 2 ) n = ±l,±2,±3 

Thus (|6]) amounts to the differential equation 
dt 



\x 



s,t) = H nss ,\x + an,s',t) (7) 



n=— 3 



Now let \ip) be a general, time-independent state vector, and express it as a 
superposition 

X 

By the requirement that be constant while the choice of basis vectors is 
changing, the time- dependence of the coefficients ip s (x,t) is determined: 

d 3 

— ip s {x,t)= H nss'^s'(x- an,t) (8) 

al n=-3 

We now evaluate the symmetries characterizing a free spin 1/2 particle. 
Translational invariance is already incorporated in t^). As for invariance 
under lattice rotations, let direction n be rotated into direction m, and let u 
be the corresponding spin rotation matrix. Then the condition for rotational 
invariance of the hopping process is 

Hmss'" = U ss'H^ s , s „U s ,} s „, (9) 

For the on-site hopping amplitudes (n = m — 0) this implies 

#cL' = e<w (io) 

For 90° rotations about the n axis, u ss i = (8 SS > +cr" s ,)/\/2 where a™ s , denotes 
a Pauli matrix. Thus, applying (^) to H^ ss , we have n = m and hence 

Hl sl = V 6 ss/ + Kal, n = l,2,3 (11) 

Cubic symmetry implies that rj and k do not depend on n. A further 180° 
rotation (u = a) about an axis perpendicular to n gives 

Hf. n)AB = V Sss> ~ k a n ss , n = 1, 2, 3 (12) 



Thus, starting from (R) we arrive at 
d 3 

— i[> a (x, t) = ei[> a (x,t) + 77^ (ip s (x + ah,t) +ip s {x - an,t)) 

dt n=l 

3 

— k a ss' (Vv + t) ~ — aw, t) ) 

n=l 

We can take the continuum limit a — > by expanding the wavefunction, 

a 2 

i/; s (x ± an, t) = ^(z, *) ± ad n ^ s (x, t) + —d 2 J) s {x, t) + C(a 3 ) (13) 
Thus we obtain 

d t ip s (x,t) = (e + 6ri)ij s (x,t)-2aKal,d n ij s ix,t)+ria 2 A4> s (x,t) + 0(a 3 ) (14) 

We now make another simple but crucial assumption — that all a priori hop- 
ping amplitudes are of the same order of magnitude, or at least not dramat- 
ically different in their order of magnitude. This means, for example, that 
hopping with a spin- flip is about as likely as without it. An argument in 
favour of this assumption is that rotational invariance on a lattice is bro- 
ken down to cubic symmetries; thus orbital angular momentum is certainly 
not conserved at the length scale of the lattice spacing. Spin is also related 
to the geometrical elements of a lattice. In the Dirac-Kahler formalism 0, 
for example, spinors arise as superpositions of 0-forms, 1-forms, etc.; in the 
lattice version of the formalism they are superpositions of sites, links, 



plaquettes, etc. 

In equation (|TJ]) our assumption means that k and 77 are of the same order 
in a. Hence, rja 2 is negligible in comparison to Ka. The terms of O (a 3 ) are 
irrelevant, too. The general solution of ( |T4]) can then be written in terms of 
Fourier components of ^ s (x, 0) as 

ip s (x,t) = e (e+6r?) * / d 3 pe ifg [cos(2a«|p1t) 5 SS > - 1 sin(2a«|plt) f> ■ a ss \ 4> s ,(p) 



Now consider the case that the Fourier amplitudes ip s (p) decay more slowly 
than exponentially with \p\ — > 00. If k had an imaginary part, there would be 
time evolution factors growing exponentially as exp(2a|Q : /t| \p\t) for some of 
the spin components. This cannot be compensated by ip s (p) by assumption, 
so the Fourier integral would diverge immediately for t > or t < 0. Hence, 
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k cannot parameterize a universal law of time evolution for all spinor wave 
functions unless it is a real number. As for the global evolution factor e*- e+6 '" t , 
an argument about superpositions of n-particle states analogous to the one 
given in the spin case shows that e + Qrj must be purely imaginary. Hence, 
that evolution factor can be absorbed by redefining the phase of the spinor 
wave function. If we finally identify the speed of light as 

c = 2aK (15) 

we obtain the free Weyl equation 

d t i/j s (x,t) = -c<j" s ,9 n ^(x,i) (16) 

One of the reasons for requiring Hamiltonians to be hermitian is to ensure 
that probabilities are conserved. For the systems considered so far, hermitic- 
ity was redundant because of the symmetries present. Let us finally consider 
a spin 1/2 particle in a situation where the symmetries are broken, but only 
weakly in next-to- leading order in a. Thus, instead of fllCf ) and ([0]), let the 
hopping amplitudes be 

Hq ss „ = e8 ss/ {8 S > S » - a^fos's"(x,t)) (17) 

H Iss" = V$ss' + (<W - a-lns's"{x, t)) (18) 

with e, rj, k of order a -1 as before and 7, the spin connection coefficients 
TT| , IT2] , |T3[ |, of order 1. If no further conditions are imposed on the 7s 



we obtain the Weyl equation in a space-time with torsion and non-metricity 
(metric-affine space-time) as currently investigated in the context of quantum 
gravity [14|]. Probabilities are still covariantly conserved if we define the 



probability current as usual [|TT], [TJ, [L3[ by 

and if we define the tensor connection coefficients T^ a by requiring the Pauli 



matrices, which are constant here, to be covariantly constant as well [11 , [12|, 



13[ . In terms of spin-matrix products this condition is 

n a a v + a» la + 7^ = 

It then follows by insertion of ([17]) and ([18]) in equation (|8|) and its complex 
conjugate that 

d*j a + r^/ = 
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The tensor metric g^ u = tr (a^ea^e) is not covariantly constant in this general 
case. Only by imposing the constraint 

sfttr(7 Q e) = a = 0,1,2,3 

with e = i(T2 the spin metric would we reproduce the Weyl equation in 



a metric space-time with torsion, expressed in an orthonormal frame JT3 
However, in view of potential applications to quantum gravity, the more 
interesting finding appears to be that the Weyl equation in a metric-affine 
space-time can be obtained on the basis of quantum-mechanical linearity 
alone, independently of scalar products for state vectors. 
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